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Abstract 

Vortices and axisymmetric vortex rings are considered in the framework of the 
subcritical nonlinear Schrodinger equations. The higher order nonlinearity present 
in such systems models many-body interactions in superfluid systems and allows one 
to study the effects of negative pressure on vortex dynamics. We find the critical 
pressure for which the straight-line vortex becomes unstable to radial expansion of 
the core. The energy of the straight-line vortices and energy, impulse and velocity 
of vortex rings are calculated. The effect of a varying pressure on the vortex core is 
studied. It is shown that under the action of the periodically varying pressure field 
a vortex ring may split into many vortex rings and the conditions for which this 
happens are elucidated. These processes are also relevant to experiments in Bose- 
Einstein condensates where the strength and the sign of two-body interactions can 
be changed via Feshbach resonance. 
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1 Introduction 

The nonlinear Schrodinger (NLS) equation 

-i^ t = l -V 2 4> + g\^, (1) 

where g = ±1, is one of the most studied equations of mathematical physics 
with applications ranging from nonlinear optics, where ip represents the elec- 
tric field, Bose-Einstein condensates (BECs) and superfluidity, where tp is the 
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field wavef unction, to fluid dynamics, where ip describes the amplitude of al- 
most monochromatic wave. The nonlinear term in this equation represents the 
action of a refraction index that depends on the electric field intensity in non- 
linear optics or corresponds to the two-body repulsive (g — —1) or attractive 
interactions (g = 1) of the bosons in superfluids and BECs. 

The formation of singularities in finite time for the focusing NLS equation 
(with g — 1) is in practice been arrested by higher order nonlinear terms 
that come from the expansion of the nonlinear refraction index in optics, from 
many-body interactions in superfluids or from the effects of the geometry of 
the trapping potential in BECs that offers the reduction in the number of 
dimensions. These leads to a modification of the focusing NLS equation to 
what is known as "subcritical" or cubic-quintic NLS equation 

-* = |vV + (siW ! -<hMV (2) 

This model was used to study cavitation and vortex nucleation (Josserand 
et al, 1995), formation of vapor droplets in superfluids (Josserand and Rica, 
1997), solitary waves in nonlinear optics (Kivshar and Agrawal, 2002), and 
bright solitons in elongated BECs (Sinha et al, 2006; Khaykovich and Mal- 
omed, 2006). 

In this paper we consider a general class of subcritical NLS (SNLS) equations 
that allow for different orders of nonlinearity 

- r0 t = ± W + ( 9l \4>\ 2 - s 2 H 2(1+7) )^, (3) 

where 7 is an integer. The choice of the parameter 7 = 1 corresponds to 
the first-order term in the expansion of nonlinearity in the correlation energy 
and results in the the cubic-quintic equation (2), whereas 7 = 3 is close to 
the value 2.8 used in density-functional theories (Dalfovo, 1992) to produce a 
quantitatively correct equation of state for superfluid helium II. 

Thus, we can view Eq. (3) to be a dimensionless form of a modified Gross- 
Pitaevskii equation 

= -|^V 2 ^ + (w^| 2{1+7) - V \ij\ 2 )ij, (4) 

where m is the mass of a boson, V is a 5— function two-body attractive inter- 
action potential and many-body repulsive interactions are characterized by a 
parameter W . Equation 4 can be written in the Hamiltonian form 

dip _ 5H 
2 



(5) 



where the Hamiltonian (energy functional) of the system is 



To find the ground state the energy functional (6) has to be minimized subject 
to the conservation of the number of particles N = J|^| 2 dx. This can be 
achieved by introducing a Lagrange multiplier \x and minimizing H — /j,N. 
The energy of the ground state becomes 

V N 2 W N 2 +~< 

and the ground state is given by t/> 2 = N/V, V = J dx. The chemical potential 
ji can be introduced into (4) explicitly by 

ip — > ■0exp(— ifit/h), (8) 

so that (4) becomes 

The ground state, ip , gives the value of the chemical potential as 

= ^o^o (1+7) " Vo^Pl (10) 



The hydrostatic pressure of the system is found as 



p = _OHo = (l + 7 )Wb 2(2+7) _ Vo 4 

dV 2 + 7 ^° 2 ^° K J 



Hydrodynamic relation for the compressibility 

1 On 



mc 2 OP 



(12) 



where n = |-?/>o| 2 is the number density, gives the expression for the speed of 
sound, c as 

c2 = (1 + i)W n 1J * - V n 
m 

We conclude that the SNLS model, in particular, can be used to study the 
effects of a negative pressure. This is not possible with the cubic NLS model 
(1) that is often used as a phenomeno logical model of superfluid helium, but 
where P ~ p 2 . 

In the last couple of decades there has been a number of experiments which 
explore the behaviour of superfluid helium at negative pressure by means of 
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ultrasound waves that produce an oscillating pressure within a small volume 
of helium (Nissen et al, 1989; Maris and Xiong, 1989; Xiong and Maris, 1989). 
In BECs experiments, Feshbach resonance is used to change the magnitude 
and sign of the scattering length. This can be modelled by Eq. (4) with a 
periodically varying Vo(t) in the presence of the external magnetic trap. In 
this paper we shall study the effects of variations of Vo on the vortex structure 
and dynamics. 



X^&X, t -> _ ,2(1+7) *' 0^000, (14) 



We cast Eq. (9) into dimensionless form by 

h 

where the healing length, b, is defined by 

h 



W mip 2 {1+1) 



(15) 



so that Eq. (9) becomes 

- ilH = ^V 2 + (H 2(1+7) - 2£|0| 2 + (1 - 20)0, (16) 
where we denoted £ = Vo/2Wo0o- 

At P = 0, n = (y (2 + 7)/2(l + 7)W ) 1/7 and the speed of sound is c 2 = (1 + 
7)7W / o n o +7 /(2 + 7)-The healing length becomes b = h/mcQ\J ^{1 + 7)/ (2 + 7). 
Since the known speed of sound in superfluid helium is approximately 238 m 



s 1 , the healing length at zero pressure is b = 0.6685y 7(1 + 7)/(2 + 7) A, 
which for 7 = 3 gives b = 1 A and the unit of time as 6.28 x 10~ 13 s. 

In dimensionless units, the wavefunction in the bulk is ip = 1, the density is 
p = \i/j\ 2 , and the local pressure and the local speed of sound are 



P = ^V +7 -£p, (17) 

c= v /(l + 7)p 1+7 -2£p. (18) 

The paper is organized as follows. In Section 2 we study the amplitude and 
energy of the straight-line vortex at various pressures (interactomic strengths). 
We find the critical value of £ at which the vortex becomes unstable to radial 
expansion for various values of 7. Section 3 is devoted to vortex rings and 
other travelling coherent structures that propagate with a fixed velocity. We 
calculate the energy and impulse of these matter waves. The time evolution of 
a line vortex and a vortex ring with a periodically varying pressure is studied 
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in Section 4. We show that the vortex ring splits into many vortex rings if 
the negative pressure reaches the critical value for the vortex instability. We 
conclude with Section 5. 



2 Vortex Lines 



A vortex line is defined by a zero of the wave function if; — 0. In cylindrical 
coordinates (r, 9, z) the wave function of the straight-line vortex takes form 

V> = /(r)exp[is0], (19) 

where s is an integer ("winding number", "topological charge"). Fluid rotates 
around the z-axis with the tangential velocity 

u = sV9 = -e g (20) 
r 



and the amplitude, /, satisfies 



l Mr d A-~f 
r dr\ dr I r 2 



+ 2£/ 3 +(l-2O/-/ 2(1+7)+1 = 0. (21) 



The boundary conditions are /(0) = and / — > 1 as r — > oo. Analysis of Eq. 
(21) shows that near the origin f(r) ~ a s r I s and at infinity 

fM „ -I * 2 * 2 (27V + 7(8 + 3, 2 ) - (8 + - 1)) 

71 J 4(l + 7 -20r 2 32(l + 7 -20 3 r 4 + ' 

(22) 

Similar to the case of the cubic defocusing NLS equation, only vortices with 
single unit of quantization, s — ±1, are dynamically stable (Josserand et 
al, 1995). Equation 21 was solved using a finite-difference discretization and 
Newton-Raphson iterations. Figure 1 gives plots of the vortex amplitudes / 
as functions of r. Figure 2 depicts the slopes at the origin of the vortex am- 
plitudes, ai, as functions of £ for various values of 7. Above a critical value of 
£crit(7) the vortex becomes unstable and expands in a radial direction as seen 
by numerically integrating Eq. (16) starting with the stable vortex and rais- 
ing the value of £. The existence of a critical pressure for the vortex stability 
has been established in the context of density-functional theory by Xiong and 
Maris (1991) (giving the value of the critical pressure as —6.5 bars) and by 
a more sophisticated theory by Dalfovo (1992) (giving the value of —8 bars). 
We observe that the instability occurs when the slope of the vortex amplitude 
at the center of the vortex becomes zero and the negative pressure forces can 
no longer be balanced by the centrifugal energy of the fluid flow. 

We can find the values of the critical pressure for the instability by evaluating 
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Fig. 1. The amplitudes of straight-line vortices for s = ±1 as solutions of Eq. (21) 
for various 7 = 1,2, 3, 4, 5 and £ = (1 + 7) /(2 + 7) that correspond to zero pressure. 
The larger values of 7 correspond to tighter vortex cores. 
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Fig. 2. The slopes at the origin, a±, of the amplitudes of straight-line vortices for 
s = ±1 as solutions of Eq. (21) for various 7 = 1,2,3,4 as functions of £. Dots 
represent numerical integration, solid lines - Eq.(23). The values of (£, ai) that 
correspond to zero pressure are shown by red crosses. 
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Eq. (17) with p 7 = Vo/2W / ^crit and calculating Wq and V at zero pressure 
using the speed of sound 238 m s" 1 , mass of the boson 6.628 x 10 -27 kg and 
the density 145.2 kg m~ 3 . For 7 = 3 we obtained the critical pressure as —6 
bar, that agrees quite well with the estimates using much more sophisticated 
density-functional theories. 

Phenomenologically, the slope at the origin obeys a simple law 

ai(0 = ^1(0)^1 -f/&rit, (23) 

where ai(0) = lim^o °i (0 an d £crit = a7 (0). Equation (23) gives the slope 
at the origin accurate to two significant digits. For instance, our calculations 
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agree with the result presented by Josserand et al (1995) done for 7 = 1. 
Numerically, we get a\ = 0.286 for £ = 5/8, and Eq. (23) gives a\ = 0.29. 
Table 1 shows the values of ai(0) and £ crit for various 7. 

Table 1. 





7 = 1 7 = 2 7 = 3 7 = 4 


ai(0) 

Ccrit 


0.9575 1.02155 1.05921 1.08409 
0.689 0.789 0.851 0.893 



To determine the energy of the vortex, we need to modify the energy functional 
(6) by the presence of the ground state. To do this we restore dimensional 
units temporarily. Following Jones and Roberts (1982) we denote by ip u the 
wavefunction of the undisturbed system of the same mass, so that 

|h/fdx = ^, (24) 

V 

where v — J v dx, and by ip^ the wavefunction of the bulk: ip — > ip^ as r — > 00. 
For the quartic term in the expression for the energy we write using Eq. (24) 

q 2 = J y H 4 rfx - = / y (|Vf - V4)^x - (^ - i>l?v. (25) 
Similarly, we obtained recursively 

q 3 = j v H 6 dx - tfv = J v m 2 - 3 rfx - - i^lfv + 3^ 2 . 

n 

-{< ~ + EM^^CW*, ( 26 ) 

k=l 

where n is a positive integer and Cf are the binomial coefficients. The terms 
(ip^e —ip^) n v are 0(\/v) and vanish as v — > 00. In this limit, the dimensionless 
energy becomes 

E 7 = i||V^| 2 rfx-£|(|^ 2 |-l) 2 rfx+^-Q 2+7 , (27) 
where Q2+7 are defined recursively as 

Q 2 = |(H 2 -l) 3 dx, 

/n 
m 2 - l) 2+ ^x+ ^(-l) fe - 1 C 2 fe +n g 2+n _ fc . (28) 
fc=l 
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Fig. 3. Values of vortex core parameter i defined by (31) of straight-line vortices 
for s = ±1 for various 7 = 1,2,3,4 as functions of £. The values of the vortex core 
parameter £ at zero pressure are —0.09 for 7 = 1, 0.22 for 7 = 2, 0.37 for 7 = 3, 
0.46 for 7 = 4. The phenomenological fit (32) is given by the solid lines. 
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In dimensional units the energy per unit length of the line vortex becomes 



E„ 



ft 2 Poo 



4:71 

2 r r°° 
+ / (R 2 -l) 2+ ^rdr 

(R 2 - l) 2+ ^- k r dr 



fdR\ 



R 2 





00 



\dr J ^ r 2 



roo 

rdr-2U (R 2 -l) 
Jo 



V dr 



7 

E 

k=l 



2+7 



(29) 



The second term in the first integral in Eq. (29) represents the classical kinetic 
energy that diverges. This can be remedied by introducing a cut-off distance 
L, corresponding to the characteristic size of the container, and writing 



00 R2 

— dr 

r 



1 i? 2 



- dr + — dr + 



r 



R 2 



dr. 



(30) 



The energy per unit length of the line vortex can, therefore, be expressed in 
the form 



E v 



^ ''Poo 
47T 



hi 



+ 



(31) 



where I can be found by numerical integration. Figure 3 shows the values of 
the vortex core parameter I as a function of £ for various values of 7. 

Phenomenologically, the relationship between £ and I has the form 

i = £crit -(T X {1- 4rit) -<T 2 (£- 4rit) 3 , (32) 
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where £ cr i t is the value of £ at £ cr i t . The values of parameters Oi and 4rit for 
various 7 are given in Table 2. 

Table 2. 

The coefficients of Eq. (32) that define the vortex core parameter I. 





7 = 1 


7 = 2 


7 = 3 


7 = 4 


^crit 


-0.299578 


0.00410 


0.15429 


0.247549 


0"1 


0.117906 


0.20718 


0.261971 


0.298006 


0-2 


0.264568 


0.449719 


0.590682 


0.701111 



3 Vortex Rings 



In this section we consider the circular vortex rings that propagate in z— direction 
preserving their form. For large vortex rings with the radius R much greater 
than the size of the core characterized by the healing length b the energy can 
be found by using the energy of the straight line vortex (31) in the region 
close to the vortex line and using only the kinetic energy term in the region 
away from the vortex line. The sum of these two contributions to energy gives 
a simple result (Amit and Gross, 1966; Roberts and Grant, 1971) 



E = -i^p^R 



ln[^]-2 + € 



(33) 



b 

The momentum of the large vortex ring becomes 

p = Kp 00 nR 2 e z , (34) 

where e z is the unit vector in the direction of the ring motion. 

The vortex rings and other travelling wave structures that propagate with a 
constant velocity correspond to stationary solutions of the SNLS equations in 
the frame of reference moving with the velocity of the ring. The wavefunction 
satisfies 

= + ( H2(1+7) ~ + 2 °)^- (35) 

We can perform a variation ip ^ ip + 5ip in the expressions for momentum 
(see Jones and Roberts, 1982) 

p = k I [{r _ 1)w - ^ ~ 1)vr] rfx (36) 
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and energy (27) and using (35) show that 5E = U5p, or 



We can differentiate Eqs. (33) and (34) with respect to R and after substitution 
into Eq. (37) obtain the expression for the velocity of the large vortex ring as 



U ' 



4nR 



, ,8R 
lul T 



(38) 



We can also use the Hamiltonian group relation (37) to derive an alternative 
form of the energy functional (27). We substitute z — > az, for a constant a in 
Eqs (27) and (36). Then using the variational relationship 



§-J(E-U P ) 



gives 

dip 



E 



I 



(39) 



2 

dx. (40) 



dz 

As the radius of the vortex ring decreases the expressions (33), (34) and (38) 
are no longer accurately describe the energy, impulse and velocity of the ring. 
The sequence of the small vortex rings and other localized disturbances can 
be found numerically by a Newton-Raphson iteration technique. The infinite 
domain is mapped by the transformations x = tan -1 (Da;) and y = tan _1 (Z)y) 
to a finite grid (— |,|)x(— |,|).Disa constant chosen to lie in the range D ~ 
0.4—0.8. The resulting equations are expressed in second-order finite-difference 
form. Taking 201 2 grid points in the finite domain, the resulting non-linear 
equations are solved by a Newton-Raphson iteration procedure using a banded 
matrix linear solver based on the bi-conjugate gradient stabilized iterative 
method. The accuracy of the obtained solutions is verified by evaluating the 
integral identities (37) and (40). 

The families of the travelling wave solutions found are qualitatively similar to 
the sequence of such waves found by Jones and Roberts (1982) for the Gross- 
Pitaevskii model. They calculated the energy E and momentum p and showed 
that the sequence in the p—E plane is given by two branches meeting at a cusp 
where p and E assume their minimum values, p m and E m . As p — > oo on each 
branch, E — > oo. On the lower branch the solutions were asymptotic to the 
large vortex rings. As E and p decrease from infinity along the lower branch, 
the solutions begin to lose their similarity to large vortex rings, and Eqs. (33), 
(34) and (38) determine E, p, and U less and less accurately. Eventually, 
for a momentum p Q slightly greater than p m , the rings lose their vorticity {ip 
loses its zero), and thereafter the solitary solutions may better be described 
as 'rarefaction waves'. The upper branch consists entirely of these and, as 
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Fig. 4. Families of travelling wave solutions of the subcritical NLS equation (16) for 
7 = 1 (red thin solid line), 7 = 2 (blue thick solid line), and 7 = 3 (black dashed 
line) and (a) £ = 1/2 and (b) £ = (1 + 7)/(2 + 7). Dots show the positions of the 
transition between a vortex ring and a rarefaction pulse on the lower branches. The 
critical velocities at which the vortex ring becomes the rarefaction pulse are given 
in Table 3. 
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p — > 00 on this branch, the solutions asymptotically approach infinitesimal 
sound waves. Figure 4 shows the positions of such cusps for the SNLS equations 
(16) for two values of £: £ = 1/2 and £ = (l + 7)/(2 + 7). The critical velocities 
at which the vortex ring becomes the rarefaction pulse are given in Table 3. 

In Berloff and Roberts (1999) we discussed how a nonlocal NLS equation with 
a higher order nonlinearity brings the vortex core parameter £ and the healing 
length b into the agreement. We observe that the same harmony is achieved 
by the subcritical NLS model (16) for 7 = 3. At zero pressure that gives 6=1 
A, the same value as the one found from the sound speed in Section 1. This 
also agrees with the experiments of Rayfield and Reif (1964), where the vortex 
ring travelling at 27 cm s _1 has the energy 10 eV. 

Table 3. 
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The critical velocity of the transition between a vortex ring and a rarefaction 
pulse. 





7 = 1 7 = 2 7 = 3 


1 

2 

1+7 
2+7 


0.7 0.91 1.03 
0.45 0.6 0.7 



4 Vortex splitting 

In this section we discuss a novel mechanism of creating vorticity in super- 
fluids. Previously, the semi-classical formation of vortices has been attributed 
to one of the four basic mechanisms: due to the existence of critical veloci- 
ties (Frisch et al, 1992), due to the transverse instabilities of dispersive waves 
(Kuznetsov and Rasmussen, 1995), due to collapse of cavities (Berloff and 
Barenghi, 2004), due to energy transfer between waves (Berloff, 2004) and 
during the condensation (Berloff and Svistunov, 2002). Here we show that 
there exists another mechanism of vortex formation based on the expansion 
and collapse of the vortex core. 

A negative pressure is generated in superfluid helium II by an ultrasonic trans- 
ducer that is used to produce periodic sound pulses (Maris, 1994). This tech- 
nique has been extremely successful in detecting electron bubbles in superflu- 
ids (Ghosh and Maris, 2005). It has been shown that the tensile strength of 
superfluid helium is much less than predicted by theory and this suggests that 
the quantized vortices may play a role in the cavitation process. 

To consider the dynamics of a straight line vortex and a vortex ring in super- 
fluid under the time- varying pressure field, we integrated Eq. (16) numerically 
forward in time. We used the 4th order accurate finite difference scheme in 
space and the 4th order Runge-Kutta in time. 

We start by considering a single straight line vortex under the action of a vary- 
ing pressure field controlled by parameter £. The pressure oscillates according 
to Eq. (17) with 

Z(t) = £ + esm(Trt/2ri), (41) 

where e is a small parameter and £o corresponds to the initial zero pressure 
£o = (l+7)/(2+7). If £o+e < £crit, then the vortex never becomes unstable and 
the size of the core increases until t — r] and decreases back to the unperturbed 
value at t = 2i] in response to the applied pressure. If £ + e > £crit> then the 
vortex core continues to increase while pressure is negative. This instability 
becomes arrested by the growing positive pressure, so the vortex core returns 
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Fig. 5. Density contour plots of the vortex under the action of periodically varying 
pressure field characterized by Eq. (41) with e = 1/20, 7 = 1 and rj = 100. The 
centre of the vortex is seen as black region, intermediate densities are shown in light 
gray, larger densities in dark gray. 




to the unperturbed value at t — At]. Figure 5 shows the density contour plots 
for the dynamics of the vortex for e = 1/20, 7 = 1 and i] = 100. The field 
around the vortex remains radially symmetric and there is no vortex splitting. 
The extra energy gained is emitted by the vortex core radially as sound waves. 

To understand the action of the pressure on the vortex we calculated the 
energy of the system which in our dimensionless units is (see Eq. (31)) 



E = 7T 



ln(L/6) + £{t) 



(42) 



Figure 6 shows the plots of £{t) as a function of t. The corresponding values 
of pressure are given by dashed line (and magnified by a factor of 100 for 
77 = 100 and by a factor of 500 for 77 = 200.). The energy of the system, 
characterized by £(t), increases between pressure minimum and maximum. To 
follow the evolution of the vortex core we plot the density of the field in the 
radial direction on Fig. 7 for e = 1/20, 7 = 1 and 77 = 100. The vortex core 
grows while the pressure stays negative (until t = 2r/), after that the core 
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Fig. 6. The vortex core parameter £(t) in (42) as a function of time for the dynamics 
of the vortex under the action of periodically varying pressure field characterized 
by Eq. (41) with e = 1/20 and 7 = 1 (solid lines). Dashed line gives the evolution of 
pressure given by Eq.(17) magnified by a factor of 100 (top panel) or 500 (bottom 
panel). 
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stabilizes and then decreases in size as pressure increases; see the top panel 
of Fig. 7. For positive pressures the core continues to decrease in size while 
having an energy density much larger than the energy of the stationary vortex 
at a given positive pressure; a large density accumulates at the vortex core at 
t = 340; see the bottom panel of Fig. 7. This extra energy is emitted by the 
vortex core as outgoing sound wave packet as it tries to regain its stationary 
density profile; see t = 350 density profile at the bottom panel of Fig. 7. The 
increase in the period of the pressure oscillations leads to qualitatively similar 
scenario. The difference is that the energy of the vortex has a longer growth 
period and, therefore, reaches larger values; see the bottom panel of Fig. 6. 

For the vortices of large radii the effect of the pressure varying field is similar 
to that on the straight-line vortex, but the velocity field created by the vortex 
loop changes the pressure distribution around the core. The vortex core col- 
lapses faster where the velocity is lower leading to the vortex splitting as Figs. 
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Fig. 7. Time snapshots of the density of a straight-line vortex under the action 
of varying pressure field characterised by Eq. (41) with e = 1/20. On the top 
panel the times are 0, 50, 100, 150, 200, 250, 300; on the bottom panel the times are 
300,310,320,330,340,350. The thinner lines correspond to later times. 
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8 and 9 illustrate for 7 = 1 and 7 = 3 respectively. In both cases we start with 
a vortex ring of radius 50 at zero pressure. The parameter £ varies according 
to Eq. (41) with e = 1/20 (for 7 = 1), e = 1/10 (for 7 = 3) and 77 = 100. 
The vortex ring splits into three rings (two of positive circulation and one ring 
with the negative circulation) on Fig. 8. There are many more vortex rings 
generated on Fig. 9. We can roughly estimate the number of vortex rings that 
will be created by calculating the available energy just before splitting. Under 
the action of the pressure field characterized by (41) with e = 1/20, 7 = 1 
and f] = 100, the single straight-line vortex emits sound when i « 8 (seen as a 
sharp cusp on Fig. (6)). The total energy available can be estimated from Eq. 
(27) and for R = 50 is about 3 times the energy of one vortex ring of radius 
50 at zero pressure. 
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Fig. 8. Density contour plots of the cross section of the vortex ring under the action 
of periodically varying pressure field characterized by (41) with e = 1/20, rj = 100 
and 7 = 1. The cross-sections of vortex rings are seen as black dots, intermediate 
densities are shown in light gray, larger densities in dark gray. An initial vortex 
of radius 506 is split into three vortex rings at the end of one period of pressure 
oscillations. 




5 Conclusions 

In summary, we considered the families of straight-line vortices and axisym- 
metric vortex rings in the framework of subcritical NLS equations with a 
focusing cubic nonlinearity and defocusing higher order term, characterized 
by the parameter 7; see Eq. (16). These models can be used to study the ef- 
fects of negative pressure in superfluids and capture the dynamics of the many 
body systems with attractive two-body interactions such as spin-polarized Li? 
(Bradley et al, 1995; Sackett et al, 1998). The strength of the two-body interac- 
tions as well as the hydrodynamic pressure are characterized by the parameter 

c 

We determined the structure and energy of the straight-line vortices and 
showed that the vortex becomes unstable to uniform radial expansion when 
the slope of the amplitude at the vortex center becomes zero and calculated 
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Fig. 9. Density contour plots of the cross section of the vortex ring under the action 
of periodically varying pressure field characterized by Eq.(41) with e = 1/10, 7 = 3 
and 7] = 100. The cross-sections of vortex rings are seen as black dots, intermediate 
densities are shown in light gray, larger densities in dark gray. An initial vortex 
of radius 506 is split into many vortex rings at the end of one period of pressure 
oscillations. 




the critical value of £ for this instability. The corresponding critical pressure 
for 7 = 3 is close to the one obtained by much more sophisticated density- 
functional theories. We calculated the energy, impulse and velocity of the 
axisymmetric vortex rings and showed the families of the travelling coherent 
structures on the impulse-energy plane. 

We considered the periodic oscillations of the parameter £ in time and eluci- 
dated its effect on the vortex core. There are three possible scenarios of the 
dynamics of the vortex. 

1. If the value of £ never reaches the critical value for the vortex instability, 
then the vortex core simply adjusts to the varying pressure. As pressure 
decreases (increases), the vortex size grows (shrinks). 

2. If the value of £ exceeds the critical value during the time evolution, then 
the vortex core grows while pressure stays negative and continues to do so 
until shortly before pressure reaches its maximum positive value. After that 
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the vortex core shrinks and the extra energy is emitted as outgoing sound 
waves. 

3. Same as in item 2, but in the presence of large variations of the velocity 
field around the core, such as in the case of the vortex rings, the vortex core 
breaks into odd number of vortices (vortex rings) to preserve the total unit 
of circulation of ±1 around the initial vortex. 

Finally, we suggest that the processes of the vortex instability and vortex 
splitting can be seen in trapped BECs where the magnitude and sign of the 
scattering length, which is represented by £ in our model, are changed by 
Feshbach resonance. In particular, our results indicate that there is a critical 
value of the interatomic attractive interactions for which a straight-line vortex 
becomes unstable. It seems likely that this instability can be arrested by the 
presence of a harmonic trapping potential. By periodically varying the sign 
and the magnitude of the two-body interactions in the trapped BEC it may 
be possible to produce the vortex spitting if the trap is not too close to being 
perfectly axisymmetric. 
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